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Abstract: The fully resummed next-to- leading-order perturbative calculation of 
the energy-energy correlation in e + e~ annihilation is extended to include the lead- 
ing non-perturbative power-behaved contributions computed using the "dispersive 
method" applied earlier to event shape variables. The correlation between a lead- 
ing (anti)quark and a gluon produces a non-perturbative 1/Q contribution, while 
non-perturbative effects in the quark- ant iquark correlation give rise to a smaller con- 
tribution In Q 2 /Q 2 . In the back-to-back region, the power-suppressed contributions 
actually decrease much more slowly, as small non-integer powers of 1/Q, as a result 
of the interplay with perturbative effects. The hypothesis of a universal low-energy 
form for the strong coupling relates the coefficients of these contributions to those 
measured for other observables. 
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1. Introduction 



The energy-energy correlation [1] (EEC) in e + e~ annihilation was one of the first 
collinear and infrared safe observables [2] for which the all-order resummation of 
perturbative (PT) radiative corrections proved to be necessary, in the back-to-back 
(as well as in the forward) kinematical configuration [3-5]. It was soon noticed that 
the comparison of the theoretical prediction with the data required the introduction 
of sizeable non-perturbative (NP) corrections. A simple model for NP effects, already 
proposed by Basham et al. in 1979 [1], suggested that they should scale as 1/Q, Q 
being the total annihilation energy (hardness scale). A more detailed model for NP 
corrections to the EEC was suggested by Collins and Soper in 1985 [6]. Operationally, 
they suggested modelling the NP effects due to the transition from partons to hadrons 
as a kind of smearing of the PT distribution. 

In recent years power-suppressed NP contributions were studied for a wide va- 
riety of hard cross sections [7]. In particular, 1/Q contributions were predicted and 
phenomenologically quantified for a number of jet shape observables such as thrust 
and jet masses, the C-parameter and jet broadenings (for a review see [8]). Follow- 
ing the technique for analysing the power-behaved contributions to hard observables 
developed in [9], we are now in a position to better understand the NP effects in 
the EEC distribution and to relate the corresponding NP parameters with those 
emerging from the analyses of jet shapes. This is the main purpose of the present 
article. 

Let us start by recalling that the energy-energy correlation (EEC) is defined as 



where the sums are over all final-state particles a and b, so that each pair of particles 
is counted twice. Here x — n ~ @a6 so that in the back-to-back region x ^ 1- 

Perturbatively, the correlation is dominated by the contribution from the primary 
qq pair, while the qg correlation produces a subleading correction. 

At the non-perturbative level, there are two physically different confinement 
effects in the EEC. The first is an additional NP contribution to the qq angular 
imbalance due to radiation of secondary gluons with transverse momenta of the 
order of Aqcd, which we call gluers [10]. This contribution scales as 1/Q 2 - A more 
important NP contribution comes from the correlation between the (anti) quark and a 
gluer with E g ~ Aqcd/Xj which scales as l/(Qx)- As a result of an interplay between 
the NP and PT effects in the EEC distribution, the naive dimensions 1/Q and 1/Q 2 
get modified, in the back-to-back region, and the power-suppressed contributions 
decrease with Q much more slowly. 





£/^EEC(x), 



(1.1) 



EEC(x) 



S(C0SX + COsQ ab ) , 
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At the parton level, the quark-quark contribution to the integrated distribution 
E(x) for small values of t — tan(^x) has the structure (to leading-logarithmic order) 

£ w -(x) = /(« s iog 2 t). 

The contribution from the (anti)quark-gluon correlation is one log down: 

Zqgix) ~ ttslogi -/(a s l0g 2 t). 

The single-logarithmic enhancement here comes from the collinear singularity of the 
qg matrix element at t = 0. However it can be effectively absorbed into the qq 
contribution. Indeed, adding the energies of the quark and the gluon(s) collinear 
with it produces the initial quark energy, so that these two terms together correspond 
to neglecting the quark energy loss in the qq correlation. Having performed the 
collinear subtraction, one is left with the residual qg contribution to S at the level of 
a correction of relative order a s . Analogously, the correlation between two secondary 
gluons starts at the a 2 level and will be neglected hereafter in the derivation of the 
resummed next-to-leading PT distribution 1 . 

As a result, the EEC at small x-, & t the perturbative level, can be simply treated 
by considering the correlation between the primary quark and the antiquark, which 
are no longer aligned, because of multiple gluon bremsstrahlung, but do not lose 
energy. 

At the NP level, the leading power-behaved contribution due to the quark-gluon 
correlation is proportional to 

s£ p) (x)oc(6)-Eg T )( X ), 

where (6), depending on Q and the angle between the two energy detectors, is the 
characteristic value of the impact parameter determining the PT distribution. In the 
back-to-back limit, x = 71 ~ & ~^ 0) one observes a power behaviour, (b) oc Q' 1 ^ 1 
with a non-integer anomalous dimension 7 [5]. As a result, the leading NP correction 
to the height of the perturbative EEC plateau at x — becomes 

dE( NP ) , dE( PT ) , ^ „ 

X = 0) oc X = 0) • . 

a cos x d cos X 

The non-integer exponent depends on the treatment of the PT coupling. In par- 
ticular, for the one-loop coupling with rif = 3(5) we obtain 7 ~ 0.32(0.36) (see 
Appendix D.3). 

The final expression for the EEC which accounts for the leading power effects 
has the following structure. After extracting a kinematical factor and the "coefficient 

lr The O (oqj PT corrections coming from the gg EEC, as well as from other sources, are taken 
care of by matching the approximate logarithmic distribution with the exact two-loop matrix ele- 
ment calculation, performed in Section 5.1. 
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function" factor C = l+O (cx B ), we are left with an expression based on the "radiator" 
lZ(b), which exponentiates in impact parameter space, 



dcosx 

T(t) = ^-Jbdb J (bQt) e~ n ^ (1 - 2bX + O (b 2 A 2 QCB )) . 



(1.2) 



The linear NP correction —2b\ originates from the quark-gluer correlation, where A 
is a parameter (with the dimension of mass) which characterises the NP interaction 
at small momentum scales, 

X = ^fM [°°dka^ P \k). (1.3) 



7T 



The issue of the PT-NP matching is explained in Sect. 5.2, and the origin of the 
Milan factor M is recalled in Appendix B. 

The radiator in (1.2) contains its own NP component which is quadratic in b, 

nib) = rc( pT >(&) + i& 2 <7, 

Cf f °° -,,2 ft Q 2 1\ rNP)/,\ (1-4) 



a = 



Strictly speaking, this contribution should have been dropped since we did not anal- 
yse a comparable quadratic effect which may come from the qg correlation. However 
we choose to keep the a effect for two reasons. Firstly, it allows us to verify that 
this quadratic NP term affects the result much less that the leading bX contribution. 
Secondly, the a contribution is logarithmically enhanced in Q, which enhancement 
should not necessarily be present in the next-to-leading power contribution from the 
qg correlation. A complete analysis of 1/Q 2 effects in the EEC remains to be done. 

In the present paper we give a comparison of available data with theoretical 
expectations based on "default" values of the relevant NP parameters, without at- 
tempting a fit to extract the optimal values. Our aim is to stimulate more detailed 
experimental studies of the EEC in the back-to-back region, where there is a partic- 
ularly interesting interplay between perturbative and non-perturbative dynamics. 



2. Kinematics and resummation 

In this section, after introducing the energy-energy correlation and the kinematics, 
we recall the relevant results of resummation in the soft limit which are needed for 
power correction studies. To this end, one needs to consider also the contributions 
coming from the reconstruction of the running coupling at large distances. They are 
obtained by using the dispersive method discussed in [9]. The specific calculations 
for the EEC are similar to the ones performed for shape variables (see [11]). They 
are described in detail in Appendix A. 
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2.1 Kinematics in the soft limit 

At the parton level, the quantity EEC receives contributions from the primary quark 
p and antiquark p and the secondary partons k^. In the soft limit the primary quark 
and antiquark belong to opposite hemispheres. Neglecting the products of energies 
Ui of the secondary partons, we have 

n 

Q 2 ■ EEC(x) = 2EE5(cos x + cos & pP ) + 4 ^ Eoji5(cos x + cos & pi ) + O (tOitOj) , 



i=i 



(2.1) 



where we have used the quark- antiquark symmetry of the matrix element. 

For the parton momenta we use the Sudakov decomposition. Introducing two 
light-light opposite vectors P and P, we write 

p = (P, p = aP + pP + p t , 1* = PiP + OiP + k ti , (2.2) 

where we have taken P along the quark direction and 2PP = Q 2 . In the soft limit 
all quantities 

a*, Pi, 1 - C, 1 - p, ^ , 

are small and of the same order, while a is much smaller (quadratic in Pt/Q)- Ne- 
glecting quadratic soft terms, we have 

EEC * <i + £>! U + c» {f ~ + 2 £(a, + A)* (> " ^ ) 



)+2x>4 2 -^)} 



(2.3) 



Q 2 P 

Here we have used oiiPi = k 2 JQ 2 , ap = p 2 /Q 2 and 

, , - 9 x , n -9 pP p t n-9 pki k t 

t = tan = tan — , tan = — — , tan ■ 



2 2' 2 Qp' 2 Qcn 

The expression (2.3) takes into account the recoil of the quark-antiquark (p ( ^ 1 
and p t 7^ 0) against (soft) secondary partons. It can be cast in a more transparent 
form by using 

P (= ^l-^ojj (i-^-S>) = !-!>-!>+••■ ( 2 - 4 ) 

where the dots correspond to terms quadratic in the soft scale. Finally, using the 
fact that the matrix element is symmetric with respect to exchange of a* and Pi, we 
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can write (apart from quadratic soft terms) 



(2.5) 



This form explicitly shows that the quantity EEC is infrared and collinear safe. In 
particular, it remains finite when a secondary gluon happens to be collinear with the 
antiquark momentum, K ti — > 0, 

-* a ■ 

Ku = hi - -jPt , (2.6) 
where the matrix element has collinear singularities (see Appendix A. 2): 

|M>( 2 W^V'4^. (2.7) 

1 1 V ipp) J 4 

Indeed, the first term on the right-hand side of (2.5) does not depend on the sec- 
ondary parton (gluon) variables. As a result, collinear and soft divergences of the 
radiation probability cancel, in the standard way, in the inclusive sum of real and 
virtual contributions. The second term is proportional to the secondary parton 
momentum, and therefore is present only in the real contribution (quark-gluon 
correlation). Here the soft singularity of the matrix element, da^ja^ is damped by 
the «i factor, while the collinear singularity, K ti — > 0, is regularised by the vanishing 
difference of the delta functions in the square brackets, the direct quark-gluon con- 
tribution to the correlation and the subtraction term due to the antiquark energy 
loss which was not included into the first term, see (2.4). 

Hereafter we shall refer to the two terms in (2.5) as the qq and qg contributions 
to the EEC, respectively. Thus reorganised, the qq contribution dominates the PT 
answer, while the qg one gives rise to the leading 1/Q NP correction. 

2.2 Resummation of soft contributions 

Resummation of multiple soft gluon radiation off the qq antenna is necessary (and 
sufficient, in the leading order) for describing the EEC in the back-to-back (small t) 
region. In this approximation the partial cross sections can be factorized as 

^ = C(a s ) dw n . (2.8) 

(7 

Here dw n stands for the normalized n soft parton emission probability, and the 
"coefficient function" C(a s ) = 1 + O (a s (Q)) is included in order to match the 
soft-resummed expression with the exact two-loop result. According to (2.5), the 
observable (1.1) acquires two contributions 

' /<T = C(a s ) ^±0! l(t) , Tit) = T qq (t) + l qg (t) , (2.9) 



a d cos x 
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where 

\ (2.10) 



^2 
Pt 

Q 2 P 2 



Zqq(t) =J2[ dWn 6 y 



(2.11) 



The distribution 1 qg {t) includes the recoiling part of the qq contribution (second term 
in the square bracket) so that, as observed before, the collinear singularities in dw n 
for k u — > are cancelled. 

In Appendix A. 2 we discuss in detail the soft parton emission probabilities dw n . 
They depend on the secondary parton momenta Ku defined in (2.6) and on the 
rescaled antiquark momentum 



p = - (2.12) 
P 



and have the form 2 



dw n = d 2 p J A^e^dW n ({^, a t }; bj , (2.13) 

where the distributions dW n are factorized in the momenta of the secondary soft 
partons. To obtain such a factorization one needs to introduce the integration over 
the impact parameter b to represent the transverse momentum conservation 

m + e ^u) = s\p+ **) = 1 70 e ^ n ^ ■ ( 2 - 14 ) 

The factorization of dW n allows the soft parton resummation. In particular one has 
£ / dW n ({k u , a,}; b) = e~^ b \ TZ(0) = , (2.15) 

with 71(b) the soft emission radiator. 

The distributions dW n are singular for n ti — > and ai — > 0. At inclusive level, 
these singularities cancel against corresponding singularities in the virtual contribu- 
tions resummed by Sudakov form factors included into dW n . As a result, the radiator 
is collinear and infrared finite. 

From (2.15) we immediately obtain the qq contribution T q q{t) (see (2.10)) 

X M = < Y j bdb UQbt) e- n{h K (2.16) 



2 the {K t i,oti\ variables are convenient for describing partons in the right hemisphere, i.e. the 
one opposite to the triggered quark, see below. 
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Notice the normalization that in the limit of no secondary emission, 1Z — > 0, one has 



The "quark-gluon" EEC, I qg (t), receives contributions from each one of the 
secondary partons (see (2.11)). Due to the factorization of dW n the sum can again 
be expressed in terms of the resummed distribution based on the standard radiator, 
with the triggered parton singled out. The details can be found in Appendix B. 

3. Soft emission radiator 

In this section we analyse the radiator, which contains both PT and NP contributions. 

The essential point is the reconstruction of the running coupling, which requires 
a two-loop analysis. To this accuracy the radiator is given by the contributions of 
one and two soft partons and has the form (see [11]) 



11(b) = / du>i(k) [1 - J (M] + / dcu 2 (hk 2 ) [1 - J (b\K tl + K t2 \)} , (3.1) 



where duj\ is the one "real" soft gluon emission distribution with one-loop virtual 
correction included; du 2 is the two non-independent "real" soft parton emission dis- 
tribution. The precise expressions for du\ and du 2 are recalled in Appendix A. 3. 
Notice that the last contribution is inclusive, i.e. the sum K n + K t2 enters as argu- 
ment of the Bessel function. 

The most natural way the running coupling appears in Minkowskian observables 
[12] is through the dispersive relation, 



where the effective coupling [9], a e f[(m), is the primitive function of the discontinuity 
of a s {m). In the PT region, m 2 ^> A 2 , the effective coupling a e s(m) differs from 
the standard a s {m) by 0(al). It is important to stress that the relation (3.2) is 
supposed to be applicable both for large and small momentum scales, and thus 
makes it possible to quantify the NP contribution to the radiator. 

By using the representation (3.2) we reconstruct the running coupling in the 
radiator and obtain the following expression, see Appendix A. 3, 




(3.2) 



Jo 



Q' 



dm 2 a cS (m) 




-d_ f Q2 dK 2 



2(1- J (bK t )]\n 




(3.3) 



J m 2 + k: 



m 2 + k 2 



First we recall the PT result and then derive the leading NP part of lZ(b). 
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3.1 PT part of the radiator 



By using (3.2) we show in Appendix A. 3 that the PT part of (3.3) reproduces the 
well known next-to- leading expression [13] 



7T 



a FT (K t )[l-J (bK t )}ln 



Q 2 e"i 



(3.4) 



Here the two-loop PT coupling a. 



PT, 



is taken in the physical "bremsstrahlung" 



scheme, in which the coupling is defined as the intensity of soft gluon radiation [14]. 
Since the observable is collinear and infrared finite, the singularity is regularized 
by the factor [1 — J (6/c t )]. In (3.4) we must keep k\ > Aq CD . 

The explicit expression for the PT radiator with the next-to-leading accuracy 
was derived in [15]. It is obtained by replacing the factor [1 — J (bK t )] by the theta- 
f unction, see Appendix C.l, 



7T 



Q 2 



. PT 



(«t) hi 



2„-o 



Q 2 e 



be™ 



(3.5) 



This gives 



+ 



16nC F 
4tt/5 V2 



— (ln(l- £)+£)- ^-ln(l-£) 



ii^(i-o + ln(1 



1 



(3.6) 



where 



P 2tt 2 



, 80 = H^_5V, ft=102-^ 



(3.7) 



and a s means a s (Q) in the bremsstrahlung scheme [14]. The first line corresponds 
to the contribution from the one-loop running coupling. The radiators with the two- 
loop and one-loop a s differ at the level of an O (af In 3 b) term which is under control 
and should be kept in the PT distributions. 

The expression (3.6) only makes sense for < £ < 1, that is, for 6 min < b < 6 max 
where 



b ■ = —e~ lE 
0mm Q e 



6 min exp 



2tt 

(3 a s 



(3.8) 



and therefore we define 



^ (PT) (6>U =oo, 

n( pT \b<b min ) = o. 



(3.9a) 
(3.9b) 
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From (2.16) we find the PT part of the qq contribution within single logarithmic 
accuracy in the soft limit, 

X^\t) ~ lP(t) = ^jhdb UQbt) e-* (PT) W . (3.10) 

We recall that the qg contribution does not contain single logarithmic PT terms. The 
matching of the approximate resummed expression (3.10) with the exact two-loop 
result will be dealt with in Section 5.1. 

3.2 NP part of the radiator 

The general expression (3.3) also contains an NP contribution. The latter is given 
by the non-analytic moments of 5a c s, the NP component of the effective coupling 
(see Sect. 5.2). According to [9], the leading NP part of the radiator, 51Z, is obtained 
from (3.3) by replacing a e s by <5a:efr an d extracting from the rest of the integrand 
the leading term non-analytic in m 2 at m 2 = 0. This term comes from the region 
k 2 ~ ttl 2 Q 2 and therefore can be obtained by expanding the Bessel function in 
(3.3), 



= b 2 ■ — - / dm 2 5a e f[(m) In 2 = \b 2 o. 

2tt J m 



m 2 + nf 

3 



(3.11) 



The upper limit in m 2 is irrelevant here since Sa e & has support at small m 2 ~ Aq CD . 
In the second line we have neglected terms which generate pieces analytic in m 2 . The 
quantity a contains two NP parameters, 

a = -A 2 ,i (hiQ 2 - |) + iA 2i2 = A 2 (lnQ 2 - i) - A! 2 . (3.12) 

Here A and A are the (log) moments of the NP effective coupling 5a C R and of its 
dispersive companion ctf p , respectively (see Sect. 5.2 below). From the relation 
between a and the NP component of the running coupling, af p , 



2tt 



£ V(ln|!-0a s NP (AO, (3-13) 



it is clear that the answer remains invariant under the choice of the scale of the 
logarithms in (3.12). 

Taking account of the NP contribution to the radiator, the full quark-quark EEC 

is 

J qq (t) = ^jbdb MQbt) e -* (PT) W . (3.14) 
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The NP effect in the quark-antiquark correlation is nothing but a Gaussian smearing 
of the PT distribution X^ PT ^(t). Indeed, introducing a two-dimensional vector t we 
can represent the answer in the form of a convolution 

T qq {t) = Q 2 / dH> G ~ {t ~ P)2%F ICPD (0 , (3.15) 
J Zira 

with X^ PT \t') the PT distribution given in (3.10). 

Equation (3.15) makes it possible to directly relate the NP parameters entering 
into the definition of a with observables describing soft hadronization. Even with 
PT radiation switched off, the direction of the leading quark undergoes a random 
walk in angle due to formation of the NP hadronic plateau. As a consequence we 
expect 

a=(k 2 ± )-n(Q), n(Q)=p h (\nQ 2 -A) (3.16) 

where (k\) is the value of the mean squared transverse momentum of primary 
hadrons in jets, ph is the density of the corresponding rapidity plateau and A the 
parameter determining the effective length of the latter. This analogy gives 

A 2 = (k 2 ± )-p h , A' 2 = A 2 {A-±). (3.17) 

A naive estimate of these numbers, ignoring the effects of resonance decays, may be 
obtained using the simplest exponential parametrization of the transverse momentum 
distribution of soft hadrons, 

P(k±) oc k±exp(-2k±/ (k±)) (k±) ~ 0.30 - 0.35 GeV , (3.18) 

together with the UA5 [16] parametrization of the charged multiplicity, 

n ch = 9.11s - 115 - 9.50 ~ 1.05 Ins- 0.39 . (3.19) 

Taking account of neutrals, we find p p ~ 1.5 and A ~ 0.4, while (Arj_) = | (k±) 2 ~ 
0.13 — 0.18 GeV 2 , so that we may expect the NP parameters to be 

^l 2 ~ 0.20 - 0.27, A' 2 ~0. (3.20) 



4. Quark-gluon correlation 

The quark-gluon distribution T qg (t) can be expressed in terms of dwi and duj 2l the 
one- and two-soft parton distributions which we have introduced for the radiator (see 
(3.1)). We have 



+ J duj 2 (kik 2 ) [11(h) + u(k 2 ) ] e »(*«+*«) I 



(4.1) 
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where, according to (2.11), the functions u(ki) which probe the EEC observable are 



The relative transverse momentum k u is defined in (2.6). The distribution dui is 
singular in the limit ot\ — > as well as when the gluon momentum becomes parallel to 
that of the radiating quark, « ti — > 0. The first (infrared) singularity is compensated 
by the a\ factor in u{k\). The collinear singularity cancels in the combination of 
delta functions in (4.2). A similar regularisation occurs in du>2 with respect to the 
"parent gluon" momentum. An additional (collinear) singularity in dui2 when the two 
offspring partons become parallel, %/«! = K^/oi2, gets absorbed into the running 
coupling determining the emission of the parent gluon, see Appendix A. 3. 

4.1 PT contribution 

As shown in Appendix C.2, the PT component of T qg (t) constitutes a small O (a s ) 
relative correction to the "quark-quark" contribution 



In the first two orders in a s this contribution is fully taken into account by merging 
the approximate resummed expression with the exact O (a^) result based on the 
matrix element calculation, as will be explained below in Section 5.1. 

4.2 NP contribution 

Hereafter we concentrate on the NP component ^ of the quark-gluon correlation, 
which is the dominant power-behaved contribution to the EEC. Notice that the soft 
approximation which has been used to derive (4.1) suffices for this purpose. 

Following the procedure introduced in [11], we compute in Appendix B the NP 
contribution and obtain 




(4.2) 



lf 9 T \t)~a s (Q).iP(t). 




(4.3) 



where the leading non-analytic piece of the trigger function is 




(4.4) 



and M the Milan factor. 



The NP contribution takes the form 




(4.5) 
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with 1qq(t) the PT distribution given in (3.10). The NP parameter A can be related 
to the first moment of the NP coupling defined Sect. 5.2 below: 

A = 2A lfi M = -A±M . (4.6) 

7T 

Using (3.10) and the relation 

the result can be expressed in terms of the mean value of the impact parameter b 
averaged over the quark distribution as follows: 

n 2 r°° 

4 NP) W = \ I bdb MQbt) ^ nb) (-26A) . (4.8) 

Eq. (4.5) has a clear physical interpretation. It describes the contribution to the 
EEC when one triggers on a gluer (a gluon with n t ~ m ~ A) in a given direction, 
t, with respect to the thrust axis. The corresponding direction of the radiating 
quark is t p — t — t g , where t g is the gluer direction with respect to the quark. This 
contribution is proportional to the gluer energy which, when expressed as the ratio 
Kt/9, produces in Eqs. (4.4) and (4.5) an extra enhancement l/t g on top of the 
standard logarithmic distribution d 2 t g /t 2 g . It is this additional singular factor which 
gives rise to the non-analytic contribution according to (4.7). 

The convolution (4.5) remains finite due to "real- virtual" cancellation. The sub- 
traction term represents the quark energy loss due to an unobserved gluer, which was 
disregarded in what we chose to call the quark-quark EEC distribution. Note that 
one consequence of this convenient subtraction convention is that what we call the 
quark-gluon contribution is not positive definite. We remark also that the structure 
of the NP quark-gluon contribution X qg does not suggest that it should be exponen- 
tiated. 

Finally, observe that in the limit in which the accompanying radiation is ne- 
glected, 71(b) — > 0, one obtains 

qg { ) Qt 3 ' adcosx Q V 2t J QsimV 1 ' 

which is the first order dispersive result, in accord with the NP expectation of [1]. 

By introducing the mean impact parameter (b) = (b) (t, Q) we can cast the NP 
qg contribution (4.8) as 

4 NP) W = -2(b)X-l qg (t). (4.10) 

For not too small values of t, such that a s log 2 t < 1, we have (b) ~ l/(tQ), which 
explains an additional 1/t enhancement of the NP term on top of the kinematical 
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Figure 1: Quasi-linear dependence log (b) (0, Q) on logQ with the expected slope 

1/t 2 factor in (4.9). In the region a s log 2 t > 1 the Sudakov suppression effects slow 
down an increase of (b) which flattens off and tends to a Q-dependent constant in 
the t — > limit. 

If we use the one-loop coupling in the (two-loop) PT radiator, this behaviour 
can be explicitly computed (see Appendix D) to yield a non-integer exponent, see 
(D.22), 

/M ,n L0894 / A Qcd V 3236 f o i- 1356 Mqcd V 3595 f 

(6)(0 ' g) " A^hrJ forn/=3 ' hrJ forn/ = 5 - 

(4.H) 

In Fig. 1 this analytical prediction for n/ = 5, shown by the dashed line, is compared 
with result of a numerical integration using the full two-loop perturbative radiator. 
The two-loop curve deviates only a little from the analytical one-loop calculation, 
which is reassuring. 

The same hadronization model that was used in the previous Section to estimate 
the parameters A2 and A' 2 gives for A± the value 

AiM = jp h (k±) ~ 0.34 - 0.40 GeV , (4.12) 

which follows from the comparison of the QCD and the "tube model" result for the 
leading power correction to the mean value of thrust, [17] 

Q{1- T) NP = 2A = 2p h (k±) ~ 1 GeV . (4.13) 
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5. Final results 



Combining the qq (2.10) and qg (4.8) contributions we obtain 

lit) = T qq - + T q9 = % rbdb.J (Qbt)e- n(PT)( V- 1 2 b2 °(l-2b\) , (5.1) 

2 Jo 

It should be clear that the NP qg correlation gives the dominant 1/Q contribution, 
while the qq effect, at the level of logQ/Q 2 , is much smaller, both formally and 
numerically. In particular, we did not consider the next-to-leading NP correction, 
potentially O (Q~ 2 ), coming from triggering qg. However, it should still be legitimate 
to keep at least the leading log Q-enhanced piece in a, provided the sub leading 1/Q 2 
correction from qg is not log-enhanced as well. To answer this question one would 
have to analyse Sfl qg further. 

5.1 Matching resummed and fixed-order predictions 

We consider the integrated EEC distribution 

w . 1 f x , da f* Atdt rfS 

EW = ^X ix T x = /„ (TTTt^ (5 ' 

and use (1.2) to derive 
S( x ) = C(a s )^ J bdb e~ n{ - b \l - 2bX) duu(l + u 2 )J Q (bQu) 



C(a s )^ J dbe- n( - b \l-2b\) 



2t 

(l+t 2 )J 1 (bQt)-—J 2 (bQt) 



(5.3) 



Neglecting corrections of the order of t 2 1 in the back-to-back region, we finally 
arrive at 

S(x) = ^(«s)y J dbJ 1 (bQt)e^ b \l-2b\) + 0(t 2 ) . (5.4) 

We now take advantage of the existing exact two-loop PT prediction for EEC. To 
this end we write 



S(x) = S rcsum ( X ) + 5S( X ) (5.5) 

where S resU m is the resummed prediction, including NP corrections, and <5£ is the 
matching correction, which takes into account additional PT contributions up to a 2 . 

In order to obtain sensible predictions at small x> we have to be careful to 
subtract and exponentiate all logarithmic terms up to this order, so that 5T, remains 
finite as x 0. The resummed expression based on the PT radiator accommodates 
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all logarithmically enhanced terms a™ log" 1 1 with m > n. The finite non-logarithmic 
correction O (a s ) is taken care of in (5.4) by the one-loop coefficient function [6] 



C(a.) = l-C F [- + -)£ 
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7T 



(5.6) 



with a s = a s (Q) in the ms renormalization scheme. At the a 2 level the first (and 
only) singular subleading logarithmic correction O (a 2 log t) appears which has not 
been taken into account by the resummation procedure. We therefore include it into 
(5.4) to define 



= C(«b) y exp 



a. 



2tt/ 



dbJ^bQt) e- n ^ b \l - 2b\) , (5.7) 



where t = tan(x/2) and r = ln(l/t 2 ). The coefficient G 2 ± was obtained by fitting the 
single-logarithmic term in the two-loop PT contribution. In numerical evaluation of 
the integral in Eq. (5.7), the condition (3.9a) was imposed, so that impact parameters 
b > &max do not contribute. We did not in fact impose the condition (3.9b) because 
its effect was found to be negligible. 
The matching correction <5£ is then 



*E(X) 



1 + (A^x) - B n r - B 12 r 2 ) 



2tt 



+ (A 2 ( X ) - B 21 r - B 22 r 2 - B 23 r 3 - B 2i r A ) (^)' 



(5.8) 



where A\ and A 2 are the one- and two-loop predictions, obtained from the program 
EVENT2 [18], and the Bi/s are the coefficients obtained by expanding Eq. (5.7) to 
second order in a s , which gives 3 



(f + vr 2 + 4C(3)) Cl + - ^ C F C A - 5 -C F n f - G 21 



B u = 


3C F 


B\ 2 = 


-c F 


B 2 \ = 


/33 


B 22 = 


(io + 


B 2 3 = 


-3C 2 


B 2 4 = 





c 2 F 



11 

— ( 

9 



71 

2 
9 1 



1) CfC * + Ts CFn > 



Requiring 8H(x) to be finite as x ~^ then gives G 2 i ^ 65. 

3 Note that terms independent of \ are irrelevant to the differential EEC and therefore we omit 
them. 
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5.2 Merging PT and NP contributions 

Within the dispersive method the analysis of the perturbative and non-perturbative 
contribution is performed by splitting the coupling into two parts 

a s (k) = a s PT (A;) + a s NP (A;) , (5.9a) 
a cS (m) = a^g(m) + 5a cS (m) . (5.9b) 

It is assumed that af p (k) has a finite support, that is, it decreases fast at large k 2 . 
This implies that a e s(m) has only non-analytic m 2 moments, 

C C dm 2 

A 2p , q = ^ / ^Sa eS (m) (m 2 y ln«m 2 , (5.10) 

with p non-integer or q ^ 0. Using (3.2) it is straightforward to relate these param- 
eters to the moments of af F (k) , 

Cf f dk 2 
2t7 



A 2p — 



(5.11) 



as follows [19] 



d q 

2p ' q ~ 



sin irp 

A 2p 

Tip 



(5.12) 



In particular one has 



A lfi = -A u A 2A = -A 2 , A 2 , 2 = -2A' 2 + 2A 2 . (5.13) 

71 

In order to define these parameters more precisely, the problem of merging the per- 
turbative and non-perturbative contributions must be addressed. The relevant pro- 
cedure was discussed in detail in [20]. It involves introducing an infrared matching 
scale pii (typically chosen to be /jj =2 GeV), above which the NP component of a s is 
assumed to be negligible. The PT prediction for a given observable contains a con- 
tribution from the region fi < fi\. If it is calculated to next-to- leading order, then the 
PT coupling is represented by its two-loop expansion with respect to a s = «ms(Q): 

^m = «. + ^(m? + £K (5.14) 

The term proportional to K accounts for mismatch between the ms and bremsstrahlung 
renormalization schemes, with K given below in (C.9). 

Defining the moments of the coupling on the interval < k < /ij, normalized in 
such a way that they would all be equal if a s (k) were constant in this region, 

a P M = H dk k p In 5 £ a s (k) , (5.15) 

q\ /ij Jo fc 
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we have 

-PT/ \ , f, Q . K q + l\ 2 

By subtraction, we can now express the non-perturbative parameters (5.11) in terms 
of the full moments (5.15). In particular we have 

M P = lh P ■ ^ [a 2p -i(^i) - «2 P T 1(^1)] , 
C 

A ' 2P = ^ ' 2^ {P lll rf [«2p-l(A*l) - «2p-l(^l)] - «2p-l,l(^l) + "2^1,1 , 

(5.17) 

where «2p-i = «2p-i,o- Note that these quantities depend, via (5.14), on the or- 
der of perturbation theory used to make the PT prediction. If this is extended to 
next-to-next-to- leading order then a further term of order a^, which can easily be 
computed, should be added to (5.14). The corresponding PT terms in (5.16), which 
diverge factorially in higher orders, represent the start of the series responsible for 
subtracting off the infrared renormalon divergence in the perturbative contribution 
to the observable. Thus there is no renormalon ambiguity in the sum of the PT and 
NP contributions. 

The NP parameters A and a introduced above are now given by 

AC 

A = M—fm [dto(fii) - a^im)] , (5.18a) 



° = {( ln /^ ~ ^) -«r T (/^i)] -«fj(/^i)| -(5.18b) 

Here M. in (5.18a) is the Milan factor resulting from the two-loop analysis discussed 
in Appendix B (see also [11]). This factor is universal for all 1/Q jet observables 
considered in e + e~ annihilation [20] and DIS processes [21] and reads 

M = l+Po 1 (2.437CU - 0.052n/) = 1.920 (1.665) for n f = 5 (0) . (5.19) 



6. Comparison with experiment 

In this section we compare the above predictions with experimental data on the 
EEC near the backward direction. At present the data are not plentiful and are 
not usually binned in the optimal way for such comparisons. Therefore, rather that 
attempting a detailed fit, we used "default" values of the relevant parameters. The 
only perturbative parameter is the QCD scale, which we fixed to be 4 

Ag| =5) = 0.23 GeV , (6.1) 

4 We consistently used ri/ = 5 in the calculation of the radiator and the matching correction ST,, 
as well as in the Milan factor where it appears more questionable. 
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which corresponds to a s (Mz) = 0.118. The three non-perturbative parameters are 
moments of the coupling a s (k) over the infrared region < k < //i, which enter into 
Eqs. (5.18) and are defined by (5.15). Choosing /ii = 2 GeV, for the first two we 
take the values 

a (2 GeV) = 0.50 , «i(2 GeV) = 0.45 , (6.2) 

which come from analyses of 1/Q effects in event shapes [8,22] and 1/Q 2 corrections 
to deep inelastic structure functions [23], respectively. For the second log-moment, 
a new parameter which has not been probed in other observables, we take the value 
according to the model of Ref. [24], which is also consistent with the values (6.2): 

ai,i(2 GeV) = 0.55 . (6.3) 

In terms of the dimensionful parameters defined in (5.11), these values correspond 
to 

AiM ~ 0.33 GeV , A 2 ~ 0.2 GeV 2 , A 2 ~ 0.0 GeV 2 (6.4) 

at Q ~ Mz- Owing to the residual Q-dependence in Eq. (5.16), A± and A 2 are 
somewhat reduced at lower energies (falling to 0.2 and 0.12 respectively at Q ~ 10 
GeV), while A' 2 remains consistent with zero. 

The theoretical predictions are compared with data on the distribution in the 
angle x — n ~ @ at a range of energies in Figs. 2 and 3. The dot-dashed curves 
show the second-order PT predictions, while the long-dashed curves display the re- 
sults of purely perturbative resummation. The short-dashed curves include the NP 
quark- ant iquark smearing effects, and the final results including the NP quark-gluon 
correlation are shown by the solid curves. 

The effect of PT resummation is to dramatically reduce the cross section at 
small X-, i- e - f° r nearly back-to-back kinematics [3-5], but not enough to match the 
data. The NP contributions give a further reduction at small x an d an enhancement 
at larger values. For Q ~ M z the NP effects are dominated by the quark-gluon 
contribution linear in b, while the quadratic NP contributions to the radiator, due 
to quark- ant iquark smearing, become important at lower energies. 

The distribution in x has a kinematical suppression of the most interesting region 
of small angles. The distribution in cosx 5 which is finite at x = 0, is more infor- 
mative. Regrettably the only data set we could find that is binned in this way is at 
the single energy Q ~ 30 GeV [25] . A comparison with the theoretical expectations, 
again using the default parameters given above, is shown in Fig. 4. 

All the predicted distributions are in reasonable agreement with the experimental 
data. We would like to stress, however, that the cos# distribution is much more 
sensitive to the NP effects. With more precise data binned in cos 9 over a wide range 
of energies, it should be possible to attempt quantitative fits to extract the values of 
the important NP parameters, including the new quantity a^i. 
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Figure 2: PLUTO [25] data on the x distribution of the EEC, compared with PT and 
NP predictions. 



7. Discussion 

In this paper we have investigated the leading power-behaved non-perturbative con- 
tributions to the EEC. In particular we have demonstrated that the power-suppressed 
contributions to the EEC distribution in the back-to-back region are strongly mod- 
ified by the interplay with purely perturbative multiparton emission effects. Thus 
the expected 1/Q behaviour of the leading non-perturbative term due to the quark- 
gluon correlation turns into 5 (1/Q) ' 32 " ' 36 , while the Q-dependence of the contribu- 
tion due to NP smearing effects in the quark-quark EEC, \ogQ/Q 2 , slows down to 
(1/Q) ' 58 " ' 65 - The latter effect should also be present in the differential transverse 
momentum distribution of massive Drell-Yan lepton pairs in hadron-hadron colli- 
sions, at small transverse momenta, p± <C M. Since the Drell-Yan process is fully 
inclusive with respect to gluons, the "1/(5" effect which was leading in the EEC 
case should be absent from the transverse momentum distribution, as it is from the 
integrated cross section [9,28]. 

5 The values of the exponents we present here correspond to n/ = 3 and 5, respectively. These 
estimates are based on an analytical treatment using the one-loop coupling. The actual two-loop 
exponent of the quark-gluon contribution, in particular, is even smaller at achievable energies, see 
Fig. 1. 
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Figure 3: SLD [26] and OPAL [27] data on the \ distribution of the EEC, compared with 
PT and NP predictions. 

At the perturbative level, the present analysis includes the fully resummed 
next-to-leading logarithmic expression for the EEC distribution based on the two- 
loop radiator, which has been matched with the exact order a"l result provided by 
EVENT2 [18]. 

As far as non-perturbative physics is concerned, the aim of this paper was to 
demonstrate consistency with the general framework provided by the dispersive ap- 
proach and with the concept of universality of confinement effects. Therefore we 
have not attempted a detailed quantitative analysis but rather have compared with 
expectations based on other processes. 

The leading NP effects are controlled by three phenomenological parameters. 
The most important of them, which determines the NP qg contribution, is the one 
that describes 1/Q contributions to the means and distributions of various jet shapes. 
The value of this parameter, a (2 GeV) ~ 0.50, we have taken from jet shape 
phenomenology. 

The other two parameters, a± and a^i, determine the log Q-enhanced and the 
constant terms of the NP "l/Q 2 " qq contribution respectively. The first of the two, 
«i(2 GeV) ~ 0.45, we have taken from the analysis of power corrections to the DIS 
structure functions. Finally, the log-moment ai^{2 GeV) ~ 0.55 we have borrowed 
from the model of Ref. [24], since it is a new quantity which has not yet been probed 
in other processes. 
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Figure 4: PLUTO [25] data on the cos# distribution of the EEC, compared with PT and 
NP predictions. 



The results of these comparisons over a broad range of energies, from 8 to 91 
GeV (Figs. 2-4), are encouraging. They show that the EEC distribution in cos 6* 
in the back-to-back region, 9 ~ 7r, is highly sensitive to NP effects, and therefore a 
fuller experimental investigation of this region would be most welcome. 

In Sects. 3.2 and 4.2 we pointed out the relation between these NP parameters, 
or rather the A's given by Eqs. (5.17), and the characteristics of the rapidity plateau 
in "soft" hadron production, Eqs. (3.17) and (4.12). The standard values of the 
mean transverse momentum, (k±) ~ 0.3 GeV, and the number density, pt — 1.5, 
give values of the NP parameters in reasonable agreement with those obtained from 
other data and from the model of Ref. [24]. 

The EEC in the back-to-back region has previously been studied theoretically 
and phenomenologically by Collins and Soper [4, 6] . As far as the perturbative 
aspects are concerned, what is new in the present paper is the complete matching 
of resummed and fixed-order predictions, including exponentiation of all logarithmic 
terms up to two-loop order. 

Concerning the non-perturbative effects, Collins and Soper were the first to point 
out the necessity of a leading NP contribution that is linear in the "impact parameter" 
b. They also estimated the coefficient of this contribution by fitting low-energy data. 
In our approach such a term arises inevitably from the quark-gluon correlation and 
its magnitude is known from other observables, in particular jet shapes. Contrary to 
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the assumption of Collins and Soper, our approach does not suggest that such linear 
terms should be exponentiated. 

The dispersive approach also gives rise naturally to contributions that are quadratic 
in b (and logQ enhanced), which can be interpreted as a NP smearing of the quark- 
antiquark correlation. Collins and Soper's parametrization allowed for such contri- 
butions but they were not included in their comparisons with experiment. In our 
treatment the linear and quadratic contributions are comparable at low energies, with 
the former becoming dominant at Q ~ Mz- This emphasises again the importance 
of comprehensive experimental studies over the widest possible range of energies. 
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A. Radiator 

A.l Phase space and momentum balance 

In terms of Sudakov variables, the phase space for the emission of the primary quark- 
antiquark pair together with n partons 

**» = (2-) 4 * 4 (p+p + X>-q) ^(p 2 )0^mU^ys^) , (A.i) 

can be written as 

Mn = d Jl c ± dht m + ^- ku) j» n ^ . (a.2) 

Introducing parton transverse momenta with respect to the antiquark direction, K ti 
(2.6), and the angular antiquark variable p defined in (2.12), we can write 

»• = % *r *V+ £ *.) CP n f 57^5 • < A - 3 > 

with 
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Since the soft matrix elements are factorized it is convenient to express also the phase 
space in a factorized form. This is obtained by introducing the impact parameter b 
to represent the transverse momentum conservation 

In the soft limit the upper bounds of the parton momentum integrations can be 
arbitrarily chosen as K ti < Q and ctj < 1. An improper treatment of the hard 
region of the phase space is then corrected by introducing the coefficient function 
factor C(a s ) (see (2.8)) and performing the matching of the approximate resummed 
expression with the exact matrix element calculation to the two-loop order, which 
was discussed in detail in Section 5.1. 

We will say that a secondary parton with 

k 2 k 2 
a i> A = — % , or A > Oii - 



a,Q 2 ' y% 1 PiQ 2 ' 

is emitted in the right- or /e/lf-hemisphere respectively. (Within this convention the 
quark belongs to the left hemisphere.) We have chosen the Sudakov representation 
based on the quark momentum direction. Therefore, as long as the quark and an- 
tiquark are generally not back-to-back, the invariant phase space is not symmetric 
with respect to left-right exchanges. 

However, (A.4) remains symmetric with respect to the R L hemispheres at the 
level of the terms linear in gluon momenta, which approximation is sufficient both 
for deriving the resummed PT distribution and for extracting the leading NP effects. 
Indeed, to this accuracy we may write 



p(= i-X> \ 1 -J2^- a ) ~ ri^-aoa-A), ( A - 5 ) 

V i=i J V i=i J i=i 

where we have neglected the quadratic terms O (onaj) in the first factor and both 
O (Pifij) and cr oc p 2 in the second factor. In conclusion, in the soft limit including 
quark recoil, we may use 



- L .- L 



da { -${ciiQ - k ti ) + d(3i — —^■&{(3 i Q - k ti ) 

Oti Pi 



(A.6) 



Here we have split the radiation into two hemispheres. For the emission in the right 
hemisphere, ctj > we included the factor (1 — «j) from (A.5), and similarly the 
factor (1 — (3i) for the emission in the left hemisphere. 
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A. 2 One-loop radiator 

The soft multi-gluon radiation probability at one loop (multiple independent soft 
gluon emission) takes the form 

dw n = d<$> n \M n \ 2 ~ d 2 p / — — 



(2tt)= 



x I TT [ Cf«s^^ c ^ 

n! 11 \ 7T vr4 



- r 1 - tti 1 - Pi 

da { -^{aiQ - k u ) + dp { — —^^(PiQ - k u 

OLi Pi 



(A.7) 



where we have expressed the phase space for the emission in the right and left hemi- 
sphere in terms of a,i and Pi respectively. Here we have substituted K t for k t in the 
theta-functions determining the Right-Left hemispheres. This pretty voluntary ac- 
tion is safe: the mismatch between the lower limits of the a-integration expressed 
in terms of the transverse momentum defined with respect the quark, kt, and the 
antiquark, K t , is relatively small for small k t and/or small p t , which is our region of 
interest . 

The factor (1 — a) /a is the classical part of gluon emission which, according to 
the celebrated Low-Barnett-Kroll theorem [29] embodies both the soft singularity, 
da/a, and the first linear correction, da ■ 0{\). Taking account of the true hard 
gluon radiation, \ada, this factor gets promoted to the full quark-gluon splitting 
function P(a), 

Izsl =, P(a) = l±OzfO!. 

a v 1 2a 

Introducing the standard subtraction to accommodate virtual contributions we arrive 
at the n-gluon emission probability, 

where dW n factorizes (for given b and p< 1) into tir and ul soft gluons emitted 
into the right- and left-hemispheres, 

dW n = dW nR ■ dW nL , 71 = fl R + n L . (A.9) 

Each of the two distributions is given by the soft factorized expression. The right- 
hemisphere distribution reads 



/ da P{a) / — y d(aQ — n t ) 

71 Jo Jo K t 

x — — dai P[ai) — — 2 Jo(6««) 

i=i ^ 



(A.IO) 
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A similar expression holds for the left-hemisphere contribution. Summing over n and 
integrating dW n one obtains 



j dW n 



e' n{b) , (A.ll) 



where 11(b) is the one-loop radiator which receives contributions from the radiation 
into both hemispheres. It is given by (see (A. 10)) 

11(b) = f da^^P(a) [ Q ^[l-J (bKt)]2-&(aQ-Kt) 

Jo 2 n Jo Kt (A. 12) 

Jo K t 71 K t 

where the factor 2 in front of the theta function accounts for the two hemispheres. We 
have 1Z(0) = which ensures the expected normalization to the total cross section, 

£ / dw n (h ...k n , b) = J e^e-W = i . (A.13) 

A. 3 Two-loop radiator 

We now consider the two-loop improvement. By using the results of [11] we can 
generalise the form of the one-loop radiator to include two-loop corrections in the 
soft region. The radiator is given by (3.1) where the one and two uncorrelated soft 
parton distributions du\ and duo^ are given by 

duo^k) = AC F daP(a) ( ^-2'd(aQ-K t )l^^ + x(nt)}, (A. 14a) 

Kt t 4tt J 

dw 2 (k 1 k 2 ) = AC F dV 2 (k u k 2 ) y M 2 (k u k 2 ) . (A.14b) 

We briefly recall here the structure of these distributions [11]. 

The first distribution duo\ describes emission of a single real gluon, with a s (0) its 
on-shell coupling and x( K t) the one-loop virtual vertex correction given by 

The dispersive representation (A.15) determines x( K t) up to a scheme-dependent 
constant, of order a s 2 . Setting this constant to zero, corresponds to the choice of 
the bremsstrahlung scheme [14]. In (A.15) we have chosen to set the ultraviolet 
integration limit at /z 2 = Q 2 , rather than /z 2 = oo, in order to have the exact inclusive 
cancellation with the real emission. The error, O (a s ), induced by such a choice is 
compensated by fixing the coefficient function appropriately. The theta function in 
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(A. 14a) selects the gluons emitted in the right-hemisphere; the accompanying factor 
of 2 takes care of the contribution from the opposite (quark) hemisphere. 

The distribution cLu 2 (kik 2 ) given in (A. 14b) corresponds to the emission of two 
soft partons with 4- momenta k\ and k 2 . The uncorrelated "Abelian" two-gluon 
emission being subtracted off, the rest can be described as the contribution from the 
radiation of a (virtual) gluon followed by its decay into qq or gg in the final state. 
The corresponding matrix element M can be found in [11]. The two-parton phase 
space dT 2 reads 

2 da d 2 n t d(p 

al 2 = dm • dz — , 

an 2n 

where a = a\ + a 2 , n t = i^n + K t2 and m 2 = (k\ + k 2 ) 2 are the "parent gluon" 
variables, while z and are the momentum fraction and relative azimuth of the two 
secondary partons, qq or gg. The following kinematical relations hold: 

— * — * 

aii = za , a 2 — (1 — z)a ; q t = ; m 2 = z(l — z)q 2 . (A. 16) 

z 1 — z 

Integrating M 2 over z and 4> one finds [11] 

/ - Sa^ft. « ~ ^hm m 2 fa - ~ *} • < *■"> 

with /3 the one-loop coefficient of the beta function. 

Performing the a-integration (over both hemispheres) of the radiator 11(b) given 
in (3.1) we arrive at 

»Q 2 a„i r>2 



11(b) =AC f J q ^ In + x(Kt) ^ [ i _ Jo{bKt) 



3 



An [ Q dm 2 ( a s ^f Ql dn 2 Q 2 e~2 ( K 2 t (K 2 t +m 2 ) \ 

+ A °i ^Ul ^ ln ^w{ 2 ^ ln ^ Po\ll-Jo(bK t )]. 

(A.18) 

The collinear divergent terms, namely the logarithmic term on the second line and 
the virtual contribution x( K t), cancel if we neglect the mismatch between the real 
and virtual phase space, which is proportional to the factor \u(k 2 /(k 2 + m 2 )). This 
mismatch vanishes as m 2 / k 2 for m 2 k 2 and, as we shall see later, proves to be 
negligible within our accuracy both for the PT and NP part of the radiator. 

We are left then with the regular contribution proportional to /3q, and we get 



[V dK 2 ^ Q 2 e~l 
Jq k 2 + m 2 k 2 + m 2 



x 



(A.19) 
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The two contributions in the curly brackets in (A. 19) combine to produce the running 
coupling. To see this we invoke the dispersive representation (3.2) for a s in terms of 
the effective coupling a e s, which gives, under the m 2 integral, 



Qs(0) 

47T 



^>-^) 2 + 0(a;U = 



An 



-d 



dm 2 



(A.20) 



This leads to the following representation for the two-loop radiator (A. 19): 

3 



1Z(b) = — - / dm 2 a e s(m) 

7T Jo 



C * f A 2 , x f 

= / dm a e g{m) \ 

n Jo Jo 



dm 2 

Q 2 



f 

Jo 



dK t r i t tu m i Q 2e 2 
[1 - J (bK t )]\n 



m 2 + k 2 



m 2 + k 2 



dn 2 t 



(m 2 + k 2 ) 2 



[\-UbK t )\ (ln-^ft-4 + 1 



m 2 + k? 



(A.21) 



Perturbative equivalence. Within single logarithmic accuracy (A.21) is pertur- 
batively equivalent to the well known expression 



C, 



7T 



dn 



n ( b ) = — -Y < K t) [ 1 - MbKt) } In 



Q 2 e 



(A.22) 



with a s taken in the bremsstrahlung scheme [14]. 

To show that (A.21) and (A.22) coincide at two-loop level, we first extend the 
m 2 -integration in (A.21) to infinity. This produces a correction of the relative order 
a s (Q) ■ k 2 /Q 2 which gives rise to a non- logarithmic correction to 1Z of the order 
0(a s (Q)), which we drop as belonging to the coefficient function. Then from the 
dispersive relation (3.2) one finds 



dm 2 a e s(m) 
( m 2 + k 2y 



In 



Q 2 e 



m 2 + kf 



1 



a s (k 2 ) ln Q 2 e~ 2 



kl 



+ 



a s (k 2 )^Q 2 e-2 



k 2 



k 2 t 



O 



dm 2 [a e g(m) — a s (k 2 )] 
k 2 t [m 2 + k 2 ) 2 

*!(k 2 )_^ 

(A.23) 



k 2 



which completes the proof. 



Irrelevance of real-virtual mismatch. One more comment is warranted con- 
cerning a mismatch between the phase space boundaries for the a-integrations of the 
real and virtual contributions leading to (A. 18). Namely, the virtual integral extends 
down to a > k, 2 , while the real one is cut off at a > k 2 + m 2 . This tiny mismatch 
is nevertheless essential for the power correction analysis. Indeed, the difference of 
a-integrals of the real and virtual contributions produces 



In 



Q 2 e 



k 2 + m 2 



-In 



Q 2 e 



ln 



K,f + m 2 



m 

2 +■■■ 



{n 2 t > m 2 ) 
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which gives rise to a logarithmically enhanced O (m 2 ) contribution coming from the 
integration region m 2 <C k 2 <C b~ 2 : 



I 



q k 2 + m 2 m 4 k 2 + m 2 2 J m 2 k 2 m 2 



dn 2 t K 2 (K 2 + m 2 ) n 2 t b 2 2 f b dn 2 K l 



Taken at face value this would undermine the analysis of the quadratic power cor- 
rection to the qq EEC, since it seems to produce a logarithmically enhanced term 
of the order of m 2 ln 2 m 2 /Q 2 . However, this non-cancellation occurs at the smallest 
kinematically allowed values of a, which correspond to the values of the complemen- 
tary Sudakov variable, j3, at the edge of phase space where /3 — 1 — O (m 2 ). From 
the phase space (A. 4) there is a suppression factor 1 — (5 (coming from (, see (A. 5)) 
which degrades a contribution potentially non- analytic in m 2 down to m 4 at least. 
Therefore we can neglect this mismatch and not to worry about the fact that in this 
region neither (A. 17) nor (A. 15), which were based on soft gluon approximation, are 
valid. 

B. Quark-gluon contributions 

By using the soft multi-parton distributions described in Appendix A, we derive here 
the quark-gluon distribution T qg (t) in the form given in (4.1). We recall that 1 qg {t) 
is obtained by integrating the partially inclusive quantity 

H = J du^h) u(h) e 1 ^ 1 + J d^ihh) [11(h) + u(k 2 ) } e^ 1+f? * 2) ,(B.la) 

u{h) = 2a t [5(t 2 - t 2 J - 5(t 2 - t 2 )] ; with t p = | , t ki = ^- . (B.lb) 

In order to make explicit the finiteness of H and reconstruct the running coupling, 
we follow [11] and introduce an "inclusive source" (probing function) u'(k) to write 



H = J dui{ki) u(ki) e^ Rn + J d^ih^u'ik) 



+ 



J du 2 (k 1 k 2 ) {[u(h) + u(k 2 ) } e fbfit - u'(k) e 1 ^'} , 



(B.2) 



where in the last two terms k — k± + k 2 , the momentum of the "parent gluon" with 
positive virtuality (mass) m. There is a freedom in the choice of the expression for 
the probing function u'(k) and transverse momentum k[ describing the contribution 
to the EEC from a massive object. The only requirement is that in the m 2 — > 
limit u'(k) exp^ibKf) should coincide with the standard massless parton contribution 

— # 

u(k) exp(ibK t ) given in (B.lb). This ensures that for a collinear and infrared safe ob- 
servable the difference vanishes for m 2 — > 0, thus making the last term in (B.2) finite 
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in the limit of collinear/soft parton splitting. The first two terms of (B.2) remain 
finite to all orders, as in the case of the radiator, due to the inclusive cancellation 
between the real and virtual corrections. 

From [11] we know that the two-loop analysis can be greatly simplified if one 
defines u'(k) by replacing the transverse momentum, k 2 , in the definition of the mass- 
less source u(k) by the transverse mass, k' 2 = n 2 t + m 2 . Following this prescription 
we define 



'(*) = 2a[5(t 2 -t /2 )-5(t 2 -t 2 p )] 



u 

l_ , v , t" -\ 

K? = K 2 + m 2 , k? = W + aj?) 2 , t> k = ^L, (R3) 

with k[ and K t parallel vectors. It is straightforward to verify that the difference of 
the probing functions in the last term of (B.2) vanishes for m 2 — > 0. 

In Appendix C.2 we show that the PT component of T qg (t) gives a next-to-next- 
to-leading contribution. Here we discuss the NP component of X qg which provides the 
dominant power correction to EEC. We remark that it is legitimate to use the soft 
approximation, upon which (4.1) is based, to analyse the 1/Q power contribution. 

In [11] it was proposed to group the terms of (B.2) into three finite contributions, 

the "naive" , the "inclusive" and the "non-inclusive" contributions we shall now dis- 
cuss. 

Naive contribution. In the first term, J , we reconstruct the running coupling as 
in the case of the radiator. This is done by combining the a s (0) term in J duj\{k)u{k) 
with the regular (3 Q part of the parent gluon source, f du2(kik 2 ) u' (k) , which emerges 
after integrating the gluon decay probability over the offspring variables, see (A. 17). 
Similarly to the case of the radiator considered above, we reconstruct the effective 
coupling by using (3.2). The power contribution is then extracted by substituting 
the NP component of the effective coupling, Sa e &, for a e &. One obtains 

zr>(*) ^ jT^Mm) (^) [ *4 H- „•) . (B.4) 



Due to our choice of k' t , the trigger function fl is a function of the transverse mass, 
k 2 + m 2 , and is given by 

n(K 2 + m 2 ) = j ^ J bdb e~ n ^J (bp) J" 2da [s(t 2 - t'l) - 5(t 2 - t 2 p )j , (B.5) 

where, in the leading (linear) approximation in m ~ k ( ~ aQ, we have omitted the 
6-dependence in the inclusive source (B.3). Invoking (3.10) we arrive at 



VL(k 2 + m 



2\ 



4 / d ~t 1 {tp) I! da [ 6 { e - ^) - 6{e - e A ' (R6) 
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where we have invented the vector t p = p/Q. Analogously one can introduce the 
vector t k with modulus t' k and arbitrary direction. Performing the integration over 
k 2 in (B.4) by parts, we obtain 6 

ir \t) = °* f°° dm 2 5a cS (m) ^ . (B.7) 
it Jq m 

As a result of the K t integration by parts we have n t = 0, and the relation (B.3), 
t' k = k[ + ap, can be expressed as 

t k = tg+t p , tg = — , (B.8) 

the direction of t g being arbitrary. 

To make explicit the m 2 dependence of the trigger function Q(m 2 ) it is convenient 
to trade the a integration for that over t g to arrive at 

m f d 2 t 

This equation has a clear physical meaning. The direction of the trigger is fixed to be 
t. We have the standard logarithmic integration, d 2 t g /t 2 , with t g the direction of the 
gluer with respect to the radiating quark, weighted by the perturbative distribution 
over the quark direction, t p — t — t g . The additional singular factor l/t g comes from 
weighting by the gluer energy proportional to the ratio m/t g , where m ~ K t ~ A is a 
typical finite transverse momentum scale determining the leading power contribution. 
We conclude with the expression 

x (NP )(t) = 2CV f dm 5acs{m) p{t) = 2^o 
(For definition of the NP parameter A lj0 see Sect. 5.2.) 

Inclusive contribution. The "inclusive" contribution T in (t) is obtained by sum- 
ming the virtual correction x in / duj\{k) u(k), see (A. 14a), together with the loga- 
rithmically divergent part of J oL^fcifa) u'(k), see (A. 17). Its NP part is 

(NP) 2C A C F f° , 2 , , J ~d\ f 00 , K 2 (K 2 + m 2 ) n 

fi in = fi(/t 2 + m 2 ) - fi(/t 2 ) . 

(B.H) 

6 Since only the terms non- analytic in m 2 give non vanishing contributions to the m 2 integral 
with 5a e ff , the actual ultraviolet limit of the n t integration, n 2 < Ktmax = O (Q 2 ) , does not matter. 
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The effective coupling here has been introduced by using the relation between a e s 
and the dispersive density of the coupling a s , cf. (3.2), 



.,. f°° dm 2 . . . , d . . n a 2 (m) 

a s (k) = - ^ p(m) , p(m) = a cS (m) = -ft — + . . . 

(B.12) 

and by performing integration by parts. 

Non-inclusive contribution. Applying to the combination J cL; 2 (&i&2) [u(ki) + 
u{k,2) — u'{k)\ in (B.2) the same procedure of replacing a 2 by the derivative of a c s, 
we present the NP part of the "non-inclusive" contribution, X ni) in the form 



KNP) 



m 



Cf /",2 c- ( ~d \ f d 2 K t d(j) , 

(t) = — J dm 5a cS (m) j J —dz- M (k u k 2 ) Q m , ^ 

Vt ni = VL(K 2 n ) + VL(k 2 2 ) - VL(k 2 + m 2 ) . 

Milan factor. It is straightforward to verify that the combinations of trigger func- 
tions which enter into the inclusive and non-inclusive contributions are proportional, 
in the linear approximation, to the same function p(t) defined by (B.9), which de- 
termines the naive contribution: 

n m = n(n 2 + m 2 ) - n(m 2 ) ~ ^ • (y^ 2 + ™ 2 - «*) , (B.i4) 
n ni = ^(4) + - q(k 2 ) ~ ^ • (« a + «« - v / ^ 2 + ^ 2 ) • (B.i5) 

Such a structure is typical for the 1/Q power corrections to various jet shapes and 
leads to the universal rescaling of the naive contribution (B.4) by the so-called Milan 
factor; for details see [20]. 

Taking account of the Milan factor, the full two-loop NP component of the qg 
contribution to the EEC becomes 

l£ p) (t)=lTHt)-M = ±p(t), \ = 2A lfi M. (B.16) 

C. Perturbative analysis of sub leading corrections 
C.l Single-log corrections to the radiator 

In order to calculate IZ^^ with next-to-leading logarithmic accuracy we introduce 
the primitive function 

*(QM = —J ih T «s (k)[ln---j , (C.l) 
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where a PT is defined here in the bremsstrahlung scheme. Writing a s = a^ T (Q), 
a = af T (/c) for brevity, we have to the required accuracy 



tn W = w(~~~ + TIT k ~) ' 

k A po \a s a Air p a J 

where 

Po = 11 -2n//3, /3i = 102 - 38n//3 . (C.3) 

Using the renormalization group equation, we can replace the /c-integration by one 
with respect to a: 

<f>(Q/k) = 167l ° F r ^ (--- + — \n as 3(3 °) 

Po Ja s <x(Po + Pia/Air) \a s a 4irp R a 8vr J 

16nC F \1 1 l a Pi f a a 1 2 a \ 3/3 a 
1 In h - — — 1 — In 1 — In — In — 

a a s a s a s 47rpo \cn s ck s 2 a s J 8% a s 



ft 

Introducing 



we have, to the required accuracy, 

a a s Arcpo 
Integrating by parts we have 

rt PT \b)= I dK t [l- J (bK t )] =^HQ/K t ) 

Jo dK t 
rbQ 

= / dxJ^^bQ/x) + 0(a s (Q)) , x = bKt, 
Jo 



(C.4) 



(C.5) 



where we have neglected the non- logarithmic correction O (a s (Q)) coming from the 
upper limit, which is taken care of by the coefficient function. Taking advantage 
of the fact that $ is a slowly varying function, we may substitute into (C.l) its 
logarithmic expansion, 

$(bQ/x) = $(bQ)-<i>(bQ)-\nx + 0(a s (b- 1 )) , f(z) = -^-f , 

a In z 

to obtain, with single-logarithmic accuracy, 

POO 

Tl^ib) = $(6Q) - $(6Q) / dxJ 1 (x)\nx + ... 

Jo 



= $(6Q) - $(6Q) (In 2 - lE ) + ... = $ ( ) + O (a B ) . 



(C.6) 
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Putting everything together, we find 



3/3 



(ln(l 



Pi 



4vr/3 V2 



1 



ln 2 fl 



ln(l 



+ 



l-£ 



8tt 



in i-e) - 



2ir 



(3 \n^ + K) (ln(l-£) + 



1 - 



where now 



t 



_ a s bQe^ E 
Po — In 



(C.7) 



(C.8) 



Here we have allowed for an arbitrary renormalization scheme and renormalization 
scale /i, so that now a s = a s (fi), in both (C.7) and (C.8). The last term in Eq. (C.7) 
takes account of the scale and scheme dependence to two-loop accuracy. In the 
bremsstrahlung scheme K = 0, whereas in the Ms scheme one has 

(67 - 3tt 2 )C a - 10n f 



K 



18 



(C.9) 



C.2 Perturbative qg contribution 

To show that the qg correlation 1 qg {t) defined in (2.11) produces, at the PT level, a 
subleading O (a s ) correction to the EEC we consider the "naive" one-gluon contri- 
bution with the running coupling reconstructed as explained above. Measuring for 
brevity all momenta in units of Q we have 



d 2 bd 2 p s P n(b) 



x a da P(a)$(a — K t ) 



C F a s 



7Y 



bdbe~ n ^ 



d 2 n t e^-^ C F a s (n t ) 

^ 2 -i)-^ 2 -? 2 ) 

dn 2 t dp 2 



k 2 — a 2 p 2 



(CIO) 



a da P(a)$(a — n t ) 



x J (bp(l - a)) Jo{bK t ) 



Since, as we shall see shortly, the n t integration is non- logarithmic, we have chosen to 
neglect the running and pulled out a s as a constant factor. Getting rid of the delta 
functions and defining the common "transverse momentum" integration variable q t 
such that q t = p for the first delta function, and q t = n t /a for the second one, we 
arrive at 



7T 

oo 



bdbe- n{h) / adaP{a) 



x 



dQt 
q 2 - t 2 



J (bq t (l - a))J {bta) - J (bt(l - a))J (bq t a) ] . (C.ll) 
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Here we have substituted oo for the actual upper limit of the q t integration, q t < 1 
(i.e. q t ^ Q), because the region q t > 1 corresponds to small impact parameters 6 < 1 
(6 < l/Q) and therefore produces a negligible contribution O (t 2 ) to the answer for 
t < 1. 

The a and g t integrals converge and produce no logarithmic enhancement. For 
the sake of simplicity we shall demonstrate this property by considering the qg con- 
tribution to the height of the "PT plateau", i.e. to the EEC distribution at t — 0. 
Setting t = we have 

X (p T)(0) = cv«s r bdbe -K( b ) f 1 adaP{a) r ^ [Jo{bqt{l _ a)) _j o{bqta)] 

71 Jo Jo Jo It 



CfO- s 

!— — — 

7T 



bdbe- n{h) / da 



■ha) ■ In 

1 I - a 



2vr 



bdbe' nib) , 
(C.12) 



where we have used antisymmetry with respect to a <-> (1 — a) to substitute 
2aP(a) = 1 + (1 - a) 2 = [2 - a(l - a)} - a => (-a) , 

and 



da a In 



a 



a 



1 

2 



We conclude that X qg amounts to a non-logarithmic O (a s ) correction to the quark- 
quark contribution (2.10): 



xr } (o) 





•- / bdbe~ n ^ = 


C F a s 


7T 


2 Jo 


7T 



(PT) 



'99 



(0) 



(C.13) 



D. Analytical estimates using one- loop coupling 

We recall that the PT part of the radiator to single logarithmic accuracy is given by 

K (pt )((j) = ^£ I" A -H a , iKt) hSd, (D.l) 

with a s the two-loop coupling in the bremsstrahlung scheme [14]. In this Appendix we 
study the EEC distribution neglecting the two- loop contribution to the running 
coupling. In so doing we lose control over single logarithmic contributions to the PT 
radiator, starting from al log 3 . On the other hand, this allows us to derive analytic 
expressions for the PT distribution and for the NP corrections to it. In particular 
we shall derive in this approximation the non-integer exponents of the Q-behaviour 
of the leading NP contributions to the EEC. 



35 



D.l Perturbative qq distribution 

The PT radiator with the one-loop a s is 



where 



(L-|)ln^ + 



(D.2) 



Lee21u^, £ = 2\n(Qb), 
and the numerical value of c is 

c = = - 16 = 0.5926 , torn; = 3 . (D.3) 
The PT evaluation only makes sense for £ < L, which implies 

b < b EE — — . (D.4) 

The exponent of the radiator (D.2) has the following Mellin representation, 

3n 



with 



3^ 



F{L) = e~2 c (^) C 5 T(l + c(L - I)) = V2^L [l + O {L' 1 ) ] . 

The theta function on the left-hand side of (D.5) is ensured by the fact that the 
integrand has no singularities in the right half-plane, so that for b > bo the //-contour 
can be moved to a — > oo. Using 

bdbMbt) (^\»l(»4f-\^, (D . 6) 



we get 



xm = ^C b db MQbt) e ~ n(PT){b) = S^i t] ' (D - 7a) 

fM = *(c+v) m [-QT) {—) • (D - 7b) 

Here we have extended the 6-integral to infinity since the integrand represented by 
the right-hand side of (D.5) vanishes for b > b . 

The inverse Mellin transform (D.7) can be formally evaluated by closing the v- 
contour around the poles of T(l — v) at v — 1 + p, with p — 0, 1, . . . For t — only 
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the pole at v — 1 contributes and we derive the PT prediction for the height of the 
EEC distribution in the back-to-back region, see [5], 

r , n x %QfF{V) /A 2 elV ln ^ e- 2 ^ / Q 2 Q 2 / A 2 ei N 
J -" (0) = 4(c+l) (^J ~ 7+T V 27rCln A^ ' A 2 

(D.8) 

Making use of this result we can rewrite (D.7b) as 

For t > all the poles {y — 1 + p with p = 0, 1 . . . ) contribute, and we obtain 



3\ cln^ 



This expansion is convergent for any (tQb ). However the series is oscillating, and 
for large (tQb ) this representation is not suitable for practical computation. 

D.2 Full qq distribution 

Taking account of the NP contribution to the radiator, the full quark-quark EEC is 
given by 

J w(*) = Y J bdb MQbt) e-* (PT) ^ e-^ = J^ +l0 ° ±-f(u,t)^t), (D.ll) 



with 



- \ 2(i/-l) 

/6 -t'V2a t-t'\ [ 2 * dcj> f°°dx 2 



_^x]T ( V{v) V 2 ^2a 
h nl \^-n)J ' X °~ t 2 • 

(D.12) 

This function has no singularities in the finite //-plane. Taking the contributions 
from the poles in f(v,i) we obtain the two equivalent expansions 

o cln C+1 +P 

i (t) -i (o) V ( boQt X P + 1 ( A2e ~ 2 ) " x (t) 

(D.13) 



cln £+1+2 

-1)* /6 gt\ 2p c+l /A 2 ef X 



) 2 V 2 ; c+l+p l Q 2 ' Fp ' 
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where 



and 



3 v c l n £±l±£±H 



(-\t>l<j) n c+l+p (A 2 ei " 



n=0 



n\ c+l+p + n\ Q 



2 



In particular, Y gives the height of the plateau at t — for the full distribution, 
relative to the PT prediction (D.8): 

v°)^(o).E ( -tr c ;t:„ (#) c "^ <*.*> 



n=0 



The exponents of successive power terms slowly increase; their magnitudes oscillate 
and decrease factorially. For rif = 3 we have numerically 



3 . 0.289 / 3 x 0.482 

,2„o\ 'A 2 e2 



(D.17) 



hg r x_ 0.307 + 0.0554 (fr)' ^ | + ... 

/ a \ 0.578 , x 9 / A \ °- 964 

-1-0»7^(A) + o, 82 (^(A) + ... 

D.3 NP gg contribution 

The NP part of T qg {t) is given by the convolution 

9 A /" cPt 9 A /" 7/ 

ff '» (t) = qJ^ *> = 52 / 2^ «") • < m8 > 

with 

r>2vr 



CW= r% /^[(l + ^ 2 + 2a;cos^ 
Jo x Jo An 



v-1 



d " 2n d<p f (l + x 2 + 2a:cos^)' y - 1 -l (1 + + 2xcos^)^ 1 _ x 



2tt 1 x 2 x 2 (^!) 

(D.19) 



The latter representation leads to the series expansion 

2 
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Observing that under the exchange k ^ u — k — 1 the first factor is symmetric and 
the second is antisymmetric, we conclude that ((u) vanishes at the positive integer 
points v — 1 + p, p — 0, 1, . . . , thus cancelling the poles of the PT function f(u, t). 
The only remaining singularities of the //-integrand in (D.18) are the poles of 

((v) at v — | + p with p — 0, 1, Evaluating the Mellin transform by closing the 

contour around these poles we get an expansion similar to that for the PT distribution 

Xqq(t), 



3 

3\ cln ^Ti- 



In particular, we immediately obtain for the leading NP correction to the plateau 
height (rif = 3) 



3 \ cln 



c+ 



51 q9 (0) _ c+l/A'eS\ + \ n _M 2 ^°- 1618 



= -0.970 — . (D.22) 



2A6 X OT -(0) c+ly Q 2 J ~" ' \Q 2 

D.4 Full EEC 

Putting together the quark-quark and the quark-gluon contributions to the EEC we 
have 

Ttotit) = l qq -(t)+l qg (t) = J ^-J(v,t) | CM + ^M)} 

-K0) f ("^W* c+1 i^K^ z (D ' 23) 



where 



c+2+p 



3% cln — r^- — , q. cln ' , . 

P c+|+p \ Q 2 j 20 c + 2 + p \ Q 2 i 



(D.24) 

Setting t = we derive the leading NP power suppressed contributions to the height 
of the plateau (rif = 3) 



/ 3 X 0-323 , 3 . 1.16 

x / A \ 0.323 / A \ 1.16 



(D.25) 



where 2^(0) is the perturbative plateau in (D.8). The first NP correction comes 
from the quark-gluon and the second from the quark-quark correlation. 
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